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By using the full density matrix approach to spectral functions within the numerical renormal- 
ization group method, we present a detailed study of the magnetic field induced splittings in the 
spin-resolved and the total spectral densities of a Kondo correlated quantum dot described by the 
single level Anderson impurity model. The universal scaling of the splittings with magnetic field 
is examined by varying the Kondo scale either by a change of local level position at a fixed tunnel 
coupling or by a change of the tunnel coupling at a fixed level position. We find that the Kondo- 
peak splitting A/T/f in the spin-resolved spectral function always scales perfectly for magnetic fields 
B < 8Tk in either of the two Tx- adjusted paths. Scaling is destroyed for fields B > IOTk- On the 
other hand, the Kondo peak splitting 5/Tk in the total spectral function does slightly deviate from 
the conventional scaling theory in whole magnetic field window along the coupling-varying path. 
Furthermore, we show the scaling analysis suitable for all field windows within the Kondo regime 
and two specific fitting scaling curves are given from which certain detailed features at low field are 
derived. In addition, the scaling dimensionless quantity A/2B and 5/2B are also studied and they 
can reach and exceed 1 in the large magnetic field region, in agreement with a recent experiment 
[T.M. Liu, et al., Phys. Rev. Lett. 103, 026803 (2009)]. 



PACS numbers: 72.15.Qm, 73.23.Hk, 73.63.Kv 

I. INTRODUCTION 

As the prototypical many-body phenomenon, the 
Kondo effect^ has attracted much attention for decades, 
and it still occupies a central role for the understanding 
of many frontier problems of condense matter physics. 
Dealing with the interaction between a localized spin and 
delocalized conduction electrons, the Kondo effect has a 
characteristic energy scal^^El which is defined as the 
Kondo temperature at which Kondo conductance has de- 
creased from its extrapolated zero-temperature height to 
half of this value. By using as dimensionless unit, 
many physical phenomena can exhibit universal scaling 
relationpHSI which offer a shortcut to grasp the intrinsic 
Kondo picture. One type of scaling analysis occurs when 
an external magnetic field is applied to magnetic impu- 
rity systems (or quantum dot (QD|^: the splitting of 
zero-bias Kondo conductance peak d epen ds universally 
only on B/T^ under a magnetic fielcP^'^. That is to 
say, in a wide magnetic field window, after rescaling treat- 
ment all splittings under different parameters 5/Tk will 
follow the same path with B/Tk as the scaling variable. 
Such magnetic field window c an exte nd in theory from 
zero field up to fields B ^ Tj<-,^ ^'^ * ^^ ^ whereas in experi- 
mental measurement J^'i^'i^ magnetic fields rarely exceed 
5Tr-. For example, a detailed investigation in Kondo 
regime, in particular for larger fields at around lOOT^, 
has been made in the work by Rosch et aL^i^ These scal- 
ing properties are believed within the Kondo regime to 
be universal and cannot be affected by system param- 
eters, although such magnetic-field-induced nonequilib- 
rium effect has not reached to the level of the equilib- 
rium Kondo effect. Very recently, Liu et at^ found large 
deviations from the conventional universal scaling analy- 
sis in the measurement of Kondo differential conductance 



peak splitting (5 of a QD device. When adjusting Tk in 
two different paths (i.e. via energy level £d or the cou- 
pling strength F), they found that the splitting 5 also 
presented two different trends with the increment of Tr- 
at a larger m agn etic field B, far beyond the traditional 
scaling theorySEHUI in Kondo effect. Such specific ex- 
perimental results indicate that the scaling characteris- 
tic may encounter breakdown in certain situations, while 
all previous theorie tP'^ '' ^^' indicate that scaling analysis is 
independent on the path of adjusting Tk even at a larger 
B. Consequently it occurs to us that: does Tr- remain 
suitable as the dimensionless scaling unit, especially in 
a larger magnetic field situation? Therefore, a careful 
theoretical examination of the scaling relation and its 
universal characteristic under an external magnetic field 
is in great need. 

In this paper, we present a detailed investigation of 
scaling in the magnetic field induced splitting of the 
Kondo resonance in a Kondo QD, whose universal char- 
acters are examined by adjusting Tk in two different 
paths (i.e.: via and F). By using the full density ma- 
trix numerical renormalization grou iPl(FDM-NRG), we 
perform precise calculations of the splitting of magnetic- 
field-induced Kondo peak in the total spectral function 
(denoted as 5) as well as in the spin-resolved spectral 
function (denoted as A). The FDM-NRG method has 
many advantages compared to the conventional way used 
in previous work^, such as the complete basis set, the 
rigorous holding of the sum rules for the spectral func- 
tion, and less sensitivity of the results to the number of 
kept states. Such new features make our results more 
accurate than those from the previous studies. 

From the numerical results, we find that the splitting 
A in the spin-resolved spectral function shows perfect 
universal scaling characteristic no matter which path is 



2 



taken to adjust Tr- as long as the magnetic field B is less 
than 8Tk- For large fields B (e.g. B > IOTr-), deviations 
to scaling are found. For the spin-averaged situation (cor- 
responding to total spectral function), small deviations of 
6 from universal scaling characteristic always exist along 
the F-varying path as long as the magnetic field is ap- 
plied. Furthermore, we also fit the scaling curves ( fg and 
/a) for both splittings 6/Tk (by averaging out to elim- 
inate small deviation) and A/Tk- The fitting curves fs 
and /a have the following characteristics: (1) The fitting 
curves fg has a threshold value at about O.STr-, while /a 
is always non-zero under a field B/Tk- (2) In the low 
field {B ~ Tk), fs exhibits the linear behavior with the 
slope coefficient 0.663 which is consistent to 2/3 from the 
Fermi liquid theory. (3) Two curves fs and /a get far 
away with one another and they do not merge into a sin- 
gle curve in the large magnetic field region. In addition, 
the dimensionless scaling quantities A/2B and 6/2B are 
also investigated and we find that they can reach and ex- 
ceed 1 when magnetic field is large enough, in agreement 
with a recent experiment.'^ 

The rest of this paper is organized as follows. In Sec. 
In] we introduce our model and the FDM-NRG method. 
In Sec. |III| we show numerical results of the magnetic- 
field-induced scaling argument. Finally, Section |TV] gives 
the conclusion. 



II. MODEL AND FDM-NRG METHOD 

We consider the single-impurity Anderson model 
(SIAM)pS, and the Hamiltonian is given by: 

B BfiQi -\- Hleads ^" H coupling 



Hleads — ^ ^kgC i^^^Cko 



H, 



coupli7ig 



(1) 



In the above Hamiltonian, the Fermion operator Ckaa de- 
notes the band states of leads with energy Ska and spin a 
(t, i) , and d(j describes the impurity states with energy 
Ed] U corresponds to the Coulomb interaction between 
two electrons with different spin embedded at the impu- 
rity site, and h.c. denotes "hermitian cojugate" . Besides, 
a indicates different leads L and i?, and t„ is the coupling 
between two subsystems. Accordingly, the hybridization 
function between subsystems is given by: 



r„ =2^^|i,|2^(c^-efc„) {a^L,R) 



F = Fi -I- F, 



(2) 



we can see that F^/^, as well as F, will be constant in 
the broadband limits (adopted henceforth). 



By using the canonical transformation: 



bkcr — { — tjlCkLcr + tLCkRa) I \J t\ + t]l (3) 

we can see only the even combination of left and right 
electron states aka couples to the local impurity state: 

Hleads = '^^kaXaO-ka + '^Skb\„bka 
ka ka 

Hcoupling = ^ y^^Xada + h.C. (4) 



ka 



where V = + tj^. 

There are many ways to solve SIAM, such as the slave- 
boson mean-field theory, the Bethe ansatz approach, the 
numerical renormalization group (NRG) method, etc, 
among which the Wilson's NRG method^^ has been 
proven to be an efficient and powerful tool to deal with 
impurity systenP^', especially to obtain its Kondo fea- 
tures. If we are only interested in the transport proper- 
ties, there are standard procedures to make use of NRG 
method: 

• discretization of continuous Hamiltonian and its 
mapping to a semi-infinite chain; 

• iterative diagonalization of the chain and yield of 
flow of many-particle levels; 

• calculation of dynamic properties, such as the spec- 
tral functions. 

The first two steps can be easily accomplished by follow- 
ing Ref. [H] , and finally the Hamiltonian ([I]) becomes: 

H = lim A-("-i)/2iJ„ 

Ho - A-1/2 I^Hdot + \/v ('^^^°" '^Ld.)^ (5a) 



(5b) 



where A is the logarithmic discretization parameter, A = 
ttY,V^S{uj - Ek) = (Tl + Tj^)/2 = F/2 is also the 

hybridization function between impurity subsystem and 
leads subsystem. Hq is the starting point of the above se- 
quence of Hamiltonians Eq. ( 5b I from which an iteration 
procedure can be established. i„, denoting the hopping 
term between two neighbor sites along the chain, has an 
exponential decreasing feature with increasing n and re- 
duces to t„ — >_|(1 + A~^)A~"/^ in the limit of large n. 

Once Eq. ([5| is obtained, there are many ways avail- 
able to derive dynamic properties by calculating impurity 
spectral function in the Lehmann representation. Within 
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the past few years, the developments of NRG method, ex- 
tending its application range to various subjects includ- 
ing the bosonicPSl and time-dependentP^l situations, are 
mainly in the subject of dynamic properties. These suc- 
cessive improvements of NRG method, including the re- 
duced density matrix (DM^, the complete set of states 
combined with the reduced density matrix idea (CFSJP^, 
and most recently the full density matrix together the 
complete set of eliminated states (FDM^i^, h ave g reat 
advantages compared to the conventional wajBESI^ gg. 
pecially the last development FDM-NRG, which is far 
ahead of the conventional method, exceeds in many as- 
pects, such as dealing with the impurity problem under 
an external magnetic field, describing spectral features at 
finite temperature with high accuracy, holding the sum 
rules of spectral function rigorously, using a complete ba- 
sis set, less sensitivity of the results to the number of kept 
states , etc. 

By using the FDM-NRG method, completeness rela- 
tion reads as followj2lJ 



1 = 



N 

E E 



1^1 )ri(^i ^1 



(6) 



where TV denotes total length of the semi-infinite chain, 
iTT-min IS the first sitc at which states are discarded. 
\k, e)^ = \k)^ {Wm+i)} (8) . . . «) {Wn)} and \l, e)^ = 
lOm <^ {km+i>} {Io-at)}, in which \k}^ and |0™ 

denote the kept and discarded states of the mth iter- 
ation shell respectively, and {|crm)} represents the set 
of 4 states in the mth site along Wilson's chain (i.e.: 
|0),| t),l t,i)); thus \k,e)^,\l,e)^ can be seen as 
kept and discarded states of the mth shell containing all 
information of whole system rather than only its first 
m sites. Therefore, such treating technique naturally in- 
volves the influence of "environment"'^^. 

The next step is to make use of complete basis of set 
Eq. ^ to solve the retarded Green's function which is 
given by: 



GABit) - -ie{t)Tr[p{A{t),B{0)}] 



(7) 



where operators A,B stand for dc,dl^ respectively, 
{A, B} = AB + BA is the standard anticommutation 
relation, and p is the full density matrix: 



EEi^" 



.-per 



n (^e| = } WnP„ 



(8) 



n le 



where £J" is the Zth (discarded) eigenstate of nth itera- 

tion step, w„ = is relative weight holding the 

sum rule ~ 1 , and p„ corresponds to the density 

n 

matrix of each single shell n: {pn)i,v = 6i,i' , in 

D 

which Zn = e is its distribution function contain- 
ing only discarded states. Noticeably, such full density 



matrix treatment, in contrast to the "single-shell approxi- 
mation" of CFS-NRG method, is the spirit of FDM-NRG 
method. 

Substituting the completeness relation Eq. Q into the 
retarded Green's function Eq. Q and following the work 
by Peter et aP^, we can easily derive final results in the 
energy space as follows: 



G 



AB 



Gab 



G\ 



AB 



G 



AB 



G 



AB 




AB 



EJ' - EJi + ir] 



[p \ii 



irj 



" i,k kd uj + EJ^ - E]' + irj 



w„Ap"nk,k 

E ^"[^ 



k',k 



.[DK] „IKD] 
^l,k • ^k'J 

uj + Ef - EJ! + ir] 

r.[DK] .[KD] 

co + E'j^- Ef + iT] 



where R 



[XY] 



X 



r|i?|s)^, in which r, s S {fc, /} de- 



note the kept {k) and discarded (I) states, superscript 
X,Y ^ {K, D} indicates 'Kept' or 'Discarded', and R 
is the fermionic impurity operator d„ or dj.. Such ma- 
trix elements can be obtained within iteration process. 
Besides, [p™'"']k,k' = J2 m{k,e\pn\k' ,e)^ are elements of 

e 

the Tith component of reduced density matrix of the shell 
m. Noticeably, Eq. ^ is an equivalent formulation of 
the FDM Green function which is identical to that in 
Ref. [E]. We rewrite this result in a similar form as in 
the work by Peter et aP^ in order to reduce the time cost 
when processing calculation; furthermore, Eq. ^ can be 
easily reduced to the result in Ref. [21] by applying the 
single-shell approximation Wn = 6n,N, thus it provides an 
intuitional comparison between FDM-NRG method and 
CFS-NRG method in the final result. 

By using Eq. we can calculate the retarded Green's 
function of impurity, and consequently all dynamic prop- 
erties can be obtained straightforwardly. For instance, 
the spectral functions can be obtained at once: 



TT " 



(10) 



such treatment yields a discrete spectral function rather 
than a smooth one, and this problem can be overcome by 
using a broadening function for each single S peak. The 
kernel function in this paper has a similar form as the 
one proposed by Weichselbaum and DelftP^ 



\/Tra\ui\ 



loglcj/cj'l 
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• exp 



J-1(UJ,U! ~ < r .log|w7wol\2l I 

j^exp|-( '^'J "' Yj , |a;'|<u;o 
S{uj-uj')^PLG-H + PG-il-H) (11) 

where a, /3, 7, are relative parameters and their values 
are fixed in practical calculation at a certain temperature 
T: a = 0.8, /3 = 2T,^ = a/4,cjo = 2T. Noticeable, since 
H{uj,u;') is independent on parameter w, thus the above 
broadening function make our spectral function hold sum 
rule identically on the algorithms, also at finite field B 
and finite temperature T. 

Since the broadening of 6 function occurs only in the 
last step of NRG method for spectral function calcula- 
tion, thus the influence of broadening on the splitting of 
Kondo-peak is actually the approximation that we use 
a smooth distribution function to replace the discrete 
d function so that we can obtain smooth spectral func- 
tions. There are many logarithmic features in Kondo ef- 
fect, thus the logarithmic Gaussian broadening function 
Plg is suitable for describing Kondo effect as well as cer- 
tain features. As the analysis in Ref. the kernel 



broadening function Eq. ( 11 ) can give the most accurate 



smooth spectral functions than ever and it also has many 
advantages such as Plg is symmetric under woo;' for 
the choice 7 — a/4. 

Noticeably, although we adopted FDM-NRG methocP^ 
here, it is equivalent to CFS-NRG methocl^ in this pa- 
per, because the investigation on magnetic-field-induced 
scaling analysis in Kondo regime shown below is totally in 
the zero temperature case where FDM-NRG method nat- 
urally becomes CFS-NRG method. Therefore, the kernel 
broadening function Eq. (11) (i.e. : mixture of Gaus- 
sian and logarithmic Gaussian) is not actually used in 
this paper, but only the logarithmic Gaussian function 
Plg works. Consequently the accuracy of our results is 
determined by Plg, and such accuracy has been enough 
to grasp Kondo picture and describe precise Kondo fea- 
tures. 

What's more, there is another technique to improve 
the accuracy: the z averaginJ^HHl^ through which we can 
remove certain oscillations in spectral functions. In this 
paper we also adopt the z-averaging treatment in order 
to increase accuracy of our results. 

Finally, by using Eqs. ([9]), (|10]) and ([n]) we can calcu- 
late the impurity spectral function A^r (uJj accurately and 
exactly, from which the Kondo physics can be derived 
with great interest. Specifically these expressions are 
very helpful for investigating the magnetic-field-induced 
scaling features shown below. 

Here we have to pinpoint one difference. In practi- 
cal experiments people usually measure how the split- 
ting of Kondo conductance peak, which is extracted in 
the conductance G versus the bias Vbias curve, varies 
with the increment of magnetic field, as done in the work 
by Liu et af^ for example. By contrast, we investigate 




FIG. 1. (Color online) (a) The spin-resolved spectral func- 
tions A-f{ijj) and j4j_(aj) and the total spectral functions A{ij) 
versus the energy uj at the magnetic field B = 3.273 x 10~^. 
(b) shows the changes of the Kondo peak positions Sf, Si, 
A|, and and the splittings 5 and A with the increment 
of magnetic field B, with the inset as its magnification in a 
low magnetic field window. The parameters of SIAM for (a) 
and (b) are the same: U = l,ed = —0.5, and F — 0.16 with 
its kondo temperature Tk ~ 1.476 x 10~^; NRG parameters: 
A = 2.5, Ns = 150, and T = 0. 



the magnetic-field-induced splitting in spectral functions 
and keep the Fermi energy of left and right leads aligned 
with each other, i.e.: Vbias = 0. The two treatments are 
equivalent to one another in studying the Kondo scal- 
ing features under a magnetic field. Furthermore, the 
splitting of Kondo peak in spectral functions is also ex- 
perimentally measurable, e.g., by using an extra weak 
probe terminaP^. 



III. MAGNETIC-FIELD-INDUCED SCALING 
ARGUMENT 

We calculate the spin-resolved spectral functions 
A^{uj), Ai{uj) and total spectral function A{uj) by us- 
ing Eqs. (pl) and plj) of FDM-NRG method. In 



equilibrium situation without an applied magnetic field, 
there are three peaks coming out in spectral functions, 
among which two peaks correspond to the energy level at 
u! = ed,£d + U, while the third one corresponds to Kondo 
resonance peak roughly at the Fermi level a; = 0. Under 
an external magnetic field B, however, the Kondo reso- 
nance peak will splits into two separate peaks as shown 
in Fig. [TJa). These peak positions, denoted as 6^,61 in 
the spin-averaged case and A^,A^ in the spin-resolved 
case, can be extracted from spectral functions at the 
maximum values which can be determined by analyzing 
derivatives. Through scanning a certain region, which 
contains two stagnation points related to the splittings, 
by using binary search method, the splitting positions 
can be localized with high precision. In our calculation, 
this precision reaches O.OOITr-. The splittings 5 and A 



are then given by S 



5t and A 



A^. In 



the particle-hole symmetry case (i.e.: Sd = —U/2), these 
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FIG. 2. (Color online) (a) and (c) represents the splittings 
A and S vs B with different Sd in the spin-resolved and total 
spectral function respectively, (b) and (d) are the results of 
(a) and (c) after scaling treatment. The splitting positions 
in the spin-resolved spectral function with different Sd are 
presented in (e), while such positions in the spin-averaged 
case are shown in (f). All insets are magnifications in a low B 
window to show certain details. All parameters are the same 
as in Fig. [T] except the energy level Ed varying from —0.5 to 
—0.25, with its Kondo temperature varying from 1.476 x 10~^ 
to 5.037 X 10"^ 

peaks are symmetry corresponding to the Fermi energy 
uj = 0. FiglTJb) shows these peak positions {S^, S^, A^, 
and A^) ana the sphttings {S and A) as the function of 
magnetic field B. From FigjIJb), we can see three impor- 
tant magnetic-field-induced features: First, the splitting 
S has an obvious threshol d va lue which has been pre- 
dicted by previous theorie^^EQl^ whereas the splitting A 
in the spin-resolved case always exists even in a very weak 
magnetic field B. Second, the S-B curve can be seen as 
linear in a low field window while the A-B curve deviates 
from this feature. Third, with increasing field B, the two 
curves do not get closer to each other, conversely they 
become far away from one another. These features will 
be investigated in detail below. 

We first find out how the scaling curves change with 
varying energy level Ed which is one path of adjusting 
Kondo scaled 

Tk = ^VfU exp[ned{ed + U)/rU] (12) 

where F is the whole width of impurity's energy level. In 



Fig. [2] (a) and (c), we can see the raw data extracted di- 
rectly from spectral functions changing gradually along 
parameter-dependent paths, either in the spin-resolved 
case or in the spin-average case. It means that if we 
change some parameters, the evolution of the splittings 
A and S vs B will go along another different path at once. 
However, after scaling treatment by using Kondo scale 
Tk as the dimensionless unit, these changes of all split- 
tings with increasing B surprisingly go along the same 
path, as seen in Fig. [2] (b) and (d). Such feature, which 
is called the scaling characteristic of Kondo effect, is the 
main subject under detailed investigation in this paper. 
From the results of Fig. 2 (b) and (d), the scahng char- 
acteristic will work well for magnetic field B < 8Tk- 
Whereas with magnetic field getting larger some devia- 
tions appear, because the system has been driven closer 
to the Kondo regime in a larger B region. When the 
magnetic field B is larger than IOTk, in which Kondo ef- 
fect has been suppressed badly and the QD system is out 
of Kondo regime, obvious deviations in both of splittings 
S and A are exhibited and the scaling theory won't be 
obeyed. 

In order to clearly show that the system indeed is out 
of Kondo regime when the magnetic field B is larger than 
IOTk, we plot the curves A^/B and A^/B versus B/Tk 
by varying F (see Fig. [S]). Here it can be seen clearly 
that the deviation (i.e.: the non- universal features) are 
pushed to much larger field B with decreasing F, and in 
addition no deviations from scaling are found in a low 
field B. As a result, the deviations from scaling found 
at larger field B > IOTk are indeed due to leaving the 
Kondo regime. In the following, we will mainly focus on 
the Kondo regime with the magnetic field B less than 
IOTk. 

Under further investigation of the scaling feature, 
two differences between the spin-resolved case and spin- 
averaged case reveal. At first the splitting S has an ob- 
vious threshold field value at around 0.5Tk, by contrast 
we can't see such threshold field in the A/Tk vs B/Tk 
curve. This is because the two separated peaks in A.f{uj) 
and A^{u}) are too close to each other that they overlap 
and can not be separated in total spectral function. The 
second difference is about the scaling feature. From (b) 
and (d) (or by comparison between the insets) we can 
see the S/Tk vs B/Tk curve presents an linear scaling 
characteristic, whereas the A/Tk vs B/Tk curve shows 
a nonlinear behavior in a low field window. 

Next we study the scaling characteristics for the posi- 
tions of the Kondo peaks A^, A^, 5^, and 6^. Fig. 2 (e) 
and (f) show these peak positions (Ao-/Tr- and Su/Tk) 
vs the magnetic field B/Tk at different energy level £d- 
In approximation, all curves can merge to one curve and 
the scaling characteristics can hold. However, in detail, 
all positions of the Kondo peak slightly shift towards 
the same direction by varying energy level, in the whole 
magnetic field region (including the weak magnetic field 
B < Tk)- The deviations are within 0.5Tk and the de- 
viations in the A^ /Tk-B /Tk curves are larger than that 
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FIG. 3. (Color online) changes of A^/_B and A^/B with the 
increment of B/Tk by varying F. The parameters are the 
same as in Fig. [l] except the coupling strength F varying 
from 0.16 to 0.24, with its Kondo temperature varying from 
1.476 X IQ-^ to 9.286 x 10"^ 



in the 5^ /Tk-B /Tk curves. In addition, due to the devi- 
ations of the Ag. and (5^. are towards the same direction, 
the sphttings A (A = A^ — A^) and 5 {5 = 5^ — 5^) can re- 
duce certain parameter-induced deviations and well keep 
the scahng feature. 

Following we adjust Kondo scale Tr- in another path 
by varying the coupling strength F, in order to check 
whether scaling characteristic also works well in this sit- 
uation. We find out some unusual details reveal in the 
inset of Fig. |4] (a) although the behavior of splittings 
with increasing magnetic field shown in Fig. H (a) looks 
much similar like the one in Fig. [2] (b) and {d). Roughly 
speaking, we can consider the scaling feature as an effec- 
tive characteristic yet. For the spin-resolved case, ignor- 
ing the larger field region where system has been driven 
closer to or out of Kondo regime, the filed-induced split- 
ting A/Tfc scales perfectly for magnetic fields all the 
time. But for the spin-averaged case, we can't ignore 
such a fact that evolution of splittings 5/Tk with the in- 
crement of B /Tk don't go along exactly the same path, 
as is seen in inset of (a). Furthermore, we find that such 
slight deviation always exists as long as the magnetic 
field is applied in system, and a larger field window is 
not necessary to reveal this feature along the F varying 
path. 

The reason for inducing such deviation can be seen 
in Fig. |4] (b) and (c) where the spectral functions un- 
der two specific coupling strength F = 0.16 and 0.24 
(denoted by using subscripts 1 and 2 respectively) are 
shown. Since the magnetic field is the only factor induc- 
ing Kondo peak splitting as mentioned above, the split- 
ting positions A^/Tr- and A^/Tr- should be dependent 
only on field B and Tk- Along the Ed varying path, we 




FIG. 4. (Color online) The variations of splittings A/Tk and 
5/Tk vs B/Tk for different coupling strength F are exhibited 
in (a), and the corresponding Kondo-peak positions A^/Tk 
and 5a /Tk in the spin-resolved case and in the spin-averaged 
case are shown in (d) and (e) respectively. In addition (b) 
and (c) present the total and spin-resolved spectral functions 
for two specific cases: F = 0.16 and F = 0.24. All parameters 
are the same as in Fig. |3] 



have seen that A/Tr- scales perfectly, and here we see 
this effect also works well along the F varying path: as 
shown in Figj4] (b) and (c), the positions /S.„/Tk of the 
Kondo peak under two different coupling strength are in- 
deed nearly the same. On the other hand, the positions 
(5^(^)/Tff in the spin-averaged case are not so luck like the 
one in spin-resolved case. Certain parameters can really 
influence specific shape of spectral function and the high 
of Kondo peak (although the Kondo-peak position keeps 
fixed) , and this makes the peak positions of total spectral 
functions shift a little towards Fermi energy (i.e.: uj = Q) 
with increasing coupling strength F as shown in Figj4](b) 
and (c). Here is in the place to pinpoint that such effect 
occurs as long as B is applied, that's to say the deviation 
from conventional scaling analysis will be observed under 
any value of external magnetic field, although this devia- 
tion is usually small (within 0.2Tx in our parameters, see 
the inset of Fig|4](a)). This result is consonant with the 
experimental work by Liu et aP^ in some aspects that 
the coupling strength can indeed induce deviation from 
scaling theory, but moreover our work indicates that such 
deviation occurs all the time along F varying path, rather 
than localized in a larger field window. 

Another possible way, which explains why F can induce 
the small deviation from scaling theory, is to investigate 
Kondo peak positions shown in Fig. |4] (d) and (e). In 
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FIG. 5. (Color online) The comparison between the fitting 
curves /a and js and numerical data under various parame- 
ters are exhibited in subfigure (a), and (e) is its magnification 
in a low field window, (b) presents the fitting curves f/\/2B 
and fs/2B and the dimensionless average data (A) /2B and 
{5)/2B, where (A) and (S) are the average values correspond- 
ing to the raw numerical data in (a) . The inset (c) is the same 
with (b) but under the logarithmic coordinate, (d) presents 
/a , fs, yA and ys while the magnetic field B is in the vicinity 
of Tk- The other parameters are the same as in Fig. [l] 



present situation, the system remains symmetry {2ed + 
U — 0) and thus these peak positions keep symmetry 
as well corresponding to the Fermi energy c^; = 0. The 
Kondo-peak positions A^/Tk in the spin-resolved case 
keep nearly the same with changing F while the magnetic 
field B < 6Tk, thus A/Tk-B/Tk curve fits the scahng 
theory well in this field window. In the spin-averaged 
case, however, the peak position S^/Tk slightly deviates 
with changing F. Now the deviation of S^/Tk and S^/Tk 
are in the oppositive direction because of the symmetrical 
peak positions. So this small deviation are not eliminated 
but enlarged in the splitting S/Tk and the scaling theory 
can not hold exactly along the F varying path. 

Since observable deviation of S/Tx vs B/Tx from scal- 
ing theory is usually small as is seen in Fig. |4] (a), 
then such small deviation can be eliminated by averaging 
out so that scaling analysis can be processed. Next we 
present a scaling analysis suitable for all field windows 
within the Kondo regime, for both the spin-resolved case 
and spin-averaged case. 

In Fig. |5] (a) we exhibit the comparison between fitting 
curves and raw data extracted directly from the spectral 
functions. Such fitting curves are given by 



fA{x) = ax'' exp(— c/x) 
fsix)=dxeM~f/x"-') + h 



with the parameters: a ~ 1.1185,6 = 1.2627, c = 
0.03992; d = 18.2913,/ = 2.8572, /i = -0.3421. From 
Fig. [5] (a) we can see the fitting curves fit numerical data 
very well in whole magnetic window, including a low field 
below the conventional threshold value Q.5Tk- Thus we 
can study some detailed features in a low field with the 
help of fitting curves. For example, the splitting 5 in 
spin-averaged case does have a threshold magnetic field 
at around 0.5Tk as is seen in (e), and this result is in 
quantitative agreement with previous theorieJ^. What's 
more, by using fitting curve /a (2:) we find out that in a 
low magnetic field the curve A/Tk vs B/Tk likes the 
power function and the first derivative of A/Tk with 
respect to B /Tk is always zero at the zero-field point. 
When magnetic field B is in the vicinity of Tr-, the curves 
/a(x) and fs{x) exhibit the linear behaviors: 



yA(a;) 



aix 



hi 

&2 



(14) 



(13) 



with the parameters: ai = 0.6632, 61 = —0.1136; 02 — 
0.7097,62 = -0.3492. In Fig. [5](d), we show the curves 
fA{x), fs{x), yAix), and ysix) together. It clearly shows 
that ys/\{x) is well consistent with fs/ Aix)- In particu- 
lar, the coefficients oi of the slope of the curve ysix) is 
0.6632, which is nearly the same as 2/3 derived through 
Fermi liquid theory.121 

On the other hand, in the high magnetic field region, 
the two curves fA{x) and fsix) do not close, whereas they 
get far away with the increment of B. This result is quite 
surprise. Generally speaking it is believed that splittings 
A/Tk and S/Tk in the spin-resolved and spin-averaged 
cases should reach the limit of Zeeman splitting under a 
large enough magnetic field, thus the two scaling curves 
fA{x) and fs{x) should get closer to one another rather 
than far away with increasing magnetic field. The key 
point for the getting far away behavior is that the peak 
position A-^ always falls in the region where its second 
derivative of spin-down spectral function with respect 
to uj/Tk is negative (see Fig. |4](b) and (c)). 

In Fig. [5] (b) (c) and (d) we sho w ot her scaling fea- 
tures investigated by previous works .^^ESl 5/25 is another 
scaling dimensionless quantity and has been predicted to 
reach the limit value 1 (i.e. the splitting of the Kondo 
peak is equal to the Zeeman splitting) in the large mag- 
netic field B limits.^ Whereas in re cent experiments 
such the limit value can be exceedecPI^. In Fig. 

(b) , we show the evolution of S/2B vs the magnetic field 
B/Tk- The numerical results clearly exhibit that 6/2B 
can reach and exc eed 1 , which is in agreement with the 
recent experiment j^ ^ ^^l In addition, in the large field 
window, the curve of S/2B-B /Tk shows the logarithmic 
behavior.iiil This behavior can clearly be seen in Fig. [5] 

(c) , in which the curve of 6/2B-ln{B /Tk) is shown. 
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IV. CONCLUSION 



In conclusion, by using the FDM-NRG method we 
study the scahng characteristics of the Kondo-peak spHt- 
tings in a quantum dot system under a magnetic field. 
Similarly as in the recent experiment two different 
ways to adjust the Kondo scale Tr-, via the energy level 
£d and the coupling strength F, are considered. Both 
splittings A and 5 of Kondo resonant peaks in the spin- 
resolved spectral function and the total spectral func- 
tion are investigated in detail. We find that the split- 
ting A/Tr- in the spin-resolved case always scales per- 
fectly for magnetic fields B < 8Tx and regardless of the 
Tif-adjusted paths. When the magnetic field B is over 
IOTr-, obvious deviations are exhibited since QD system 
has been out of the Kondo regime. On the other hand, 
S/Tx in the total spectral function, which can be related 
to the experiments by Liu et aP^, deviates slightly from 
the conventional scaling theory along the F- varying path. 
Such result is consonant with the experimental work by 
Liu et al in some aspects, i.e.: the coupling strength can 
indeed induce deviation from scaling theory, but more- 



over our work indicates that such deviation occurs all 
the time as long as an external magnetic field is applied. 
Therefore Tj^ is an unsuitable parameter as the scaling 
dimensionless unit in this situation. Since the deviation 
in the splitting 6/Tk is usually small, we can still make 
the scaling analysis for 5/Tk- The fitting curves {fs and 
/a) for both splittings 5/Tk and A/Tr- are presented. 
The fitting curves fs has a threshold value at around 
O.STr", but /a is always non-zero while under the field 
B/Tk- In the low field {B ~ Tr-), fg exhibits the linear 
behavior with its slope coefficient 0.663, consistent with 
the value of 2/3 from the Fermi liquid theory. On the 
large field side, two curves fs and /a get apart, and both 
/S./2B and 6/2B can reach and exceed 1, which is also in 
agreement with a recent experiments. 
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